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Abstract 

The influence of the phases of tunneling matrix elements on the rate of the 

elastic co-tunneling at an ultrasmall normal-conducting double-junction is 

studied in a simple quantum-hole approach at zero temperature. The results 

are compared with experimental data. 
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It is known that besides ordinary tunneling of single charges through ultrasmall double- 
junctions there is the effect of macroscopic quantum tunneling (q-MQT) via virtual inter- 
mediate states Though the current of q-MQT is much smaller than that of ordinary 
tunneling, it can be observed for voltages below the Coulomb-blockade, where ordinary tun- 
neling is impossible || . Here we consider only the so-called elastic channel of q-MQT that 
is qualified by a nonvanishing conductivity at zero voltage as predicted theoretically Averin 
and Nazarov [|TJ, f|| have presented a theory of elastic q-MQT basing on |§] that takes into 
consideration the phases of tunnel matrix elements yielding the following expression for the 
conductivity G e \ at low voltages V and temperatures T 
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In this paper we try to give an modified theory of elastic q-MQT where the central 
electrode is approximated by an ID quantum hole. This is justified because of its size. We 
study a simple double-junction, consisting of two ultrasmall tunnel junctions connected in 
series and possessing capacities Ci, C2 (Cs = C\ + C 2 ), driven by the voltage V. It is 
assumed that the external electrodes are metallic bulks. Finally, the result is a formula for 
the conductivity including the discrete island level spacing Ao at the Fermi edge. Note, that 
in the continuous case A5" 1 corresponds to the energy density of states at the Fermi level. 
Gx t 2 denote the individual tunnel conductivities of the corresponding junctions. E\ and E 2 
are the electrostatic energy differences connected with tunneling at the respective electrodes 
and read as 
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n e is the integer- valued charge on the central electrode due to tunneling. In terms of a first 
order Taylor expansion in the supplied voltage V the elastic conductivity G e i turns out to 
be at low temperatures |L|, 
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Here / denotes the Fermi function and D the bulk electron density of states. The subscript 
"0" indicates the Fermi level. In the coordinate representation the tunneling amplitudes T 
are expressed by the wave functions ip in the following way 

T km = J d 3 yd 3 zT(y,z)r k (y)ip m (z). (4) 

Note, that the subscripts k and I and the coordinates y belong to the central electrode of 
the double-junction, whereas the subscripts m, n and the coordinates z are used to describe 
the banks. 

Following JTJ, [0] the conductivities of the single-junctions in Golden-rule type estimation 
are given by 

Gl , 2 = f£|KM>|»£, (5) 

where D is the density of states of the banks at Fermi level. In coordinate representation 
this looks like 

G, = ^fj d 3 y 1)2 d 3 z 1)2 T( 1 )(y 1 ,zO^ (1> (y2^2)^o(2/i)^o(^i)^o(^)V'o(y 2 )^, (6) 

and respectively in the case of the other junction. 

Owing to the fact that the tunneling matrix elements T(y, z) drop exponentially outside 
of the vicinity of the junctions @, the integrations with respect to y contribute only for 
fixed values 

2/1,2 = — g 2/3,4 = -• (7) 

Then, G e \ simplifies to 
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Using real wave functions according to |2| G e \ may be factorized in this approximation with 
regard to k and k', resulting in 
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Now, the estimation of (^j is done within a quantum hole approximation. This is different 
to the original approach of Averin and Nazarov. The central electrode is taken as ID 
quantum hole with the length L and the energetic depth E (counted from the Fermi level, 
E > 0). 

Now, the solution of the quantum mechanical problem is done within the variables [0] 
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In these terms the wave functions are 
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The finite discrete energy spectrum of the quantum hole results from the solutions of the 
transcendental equations 
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for symmetric (+) and antisymmetric (— ) wave functions, respectively. Due to the low 
conductivities of the single junctions the wave function of the island is well approximated 
by the wave function of the quantum hole. 

According to 0, ohmic behaviour of the current can be expected in case of eV <C E, 
only. Hence, we simplify the model in terms of k K which corresponds to e V <C E via 
an uncertainty relation. In this approximation the solutions of Eq. (|TT|) and fll2 ) read as 
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with tin ~ 2m*^2m*(E — Ep)/Ti. n E denotes the total number of states in the quantum 
hole of the energetic depth E. Furthermore, the following Taylor expansion in 1/(kL) of the 
required wave function values can be obtained : 
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Fig. |l| shows a comparison of different order contributions to ip(L/2)tfj(— L/2) at e(k n ) = Ef- 
By means of (0) and the approximation fli~3|) Eq. @ is given for zero temperature by 
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Here, labels the number of states below the Fermi level {np < ^e). The contributions 
of neighboring states possess different signs due to the different symmetry of the respective 
wave functions. This effect is known from [0]. Formula ( |i~5"D is our main result. It is much 
more complicated than expression (||) because the discrete energy spectrum of the central 
electrode (quantum hole) has been taken into account. 

In discussion we are going to compare our results with the experimental data of [Q. For 
the evaluation numerical values of np and Ep are necessary. According to the experimental 
setup, the values are chosen as np = 5.9 x 10 22 cm~ 3 and Bp = 5.51eV 0. The quasi-particle 
mass m* is approximated by the free electron mass m* « m ~ 9.1 x I0~ 31 kg. The parameter 
L corresponding to the diameter of the central electrodes is given by the experiment ||. The 
second parameter E describing the barrier height of the quantum hole has to be determined 
in comparison with the experimental data. Fig. |2| and Fig. |3] show the computed elastic 



conductivity G e i of sample A and B of in dependence on E, respectively. The observed 
steps arise from the occupation of a new state in the quantum hole. The parameter E is 
given by the figures where the computed .E-dependent conductivity G e \ coincides with the 
experimental value (independent of E). Due to the fluctuations of the numerical curves a 
least-squares fit was applied to them. This yields the corresponding parameters E of the 
samples A and B; Ea ~ Eb ~ 2 ... 3 eV (cf. Fig. |2] and Fig. ||]). This is in good agreement 
with earlier estimations ||. 

The pros of this model are the simple theory behind it. In addition to [g] our result 
depends on properties of the island characterized by E and L. According to the applied 
model, A ~ 0.386^, 1.12eV is found for sample A and B, respectively. These values exceed 
drastically the data concluded from the experiment in terms of the theory 0] (A = 
0.0A5eV, 0.126ey). This is not surprising because they are model-dependent parameters. 

Our approach is based on the very simple model of an ID quantum hole, which is far from 
being general due to the neglect of two dimensions and the idealization of the energetic shape 
of the barriers. Regarding to the second point it is unlikely within eV -C E that the shape 
of the upper end of barriers influences the current dramatically. An influence on the steps 
observed in Fig. ^| and Fig. [3| could be expected. This approach could be improved by using 
more sophisticated models of the barrier, as presented for instance by f|, [f|. Additionally, 
instead of the bulk values for Ep and np as given above more appropriated values should be 
taken from the experiment. 
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FIGURES 

FIG. 1. The evaluation of the different orders of a Taylor expansion of the wave function in 
a deep quantum hole. - The solid line corresponds to the numerical solution whereas the other 
curves display the Taylor expansion of given order. 

FIG. 2. The comparison of our formulae with the experimental data of sample A. - The figure 
shows the dependence of the elastic conductivity G e i on the barrier height E. We compare the 
numerical solution (solid line) with the fitted value out of the experimental paper (horizontal line). 
The steps indicate the occupation of a new state within the quantum hole. The straight line 
corresponds to a G e \ = aE b least-square fit of the numerical data. The parameters of the sample 
are L = 3nm, C\ = 0.57 x 1(T 18 F, C 2 = 3.4 x 1(T 18 F, Gi = l/(3.5Mfi), and G 2 = l/(21Mfi). 

FIG. 3. The corresponding comparison of data like Fig. 1 for sample B. Its parameters are 
L = 0.8 nm, C\ = 0.095 x 10~ 18 F, C 2 = 0.066 x 10~ 18 F, Gi = 1/(0.42AK2), and G 2 = l/(630MfT). 
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